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TECHNICAL MEMORANDUM X-64661 


BOUNDARY LAYER LOSS SENSITIVITY STUDY USING 
A MODIFIED ICRPG TURBULENT BOUNDARY 
LAYER COMPUTER PROGRAM 

SUMMARY 


The original JANNAF Turbulent Boundary Layer computer program 
was modified to use accurate boundary layer edge conditions in the calculation 
process instead of the internally calculated properties based upon ideal gas 
relationships. The Prandtl number, assumed to be constant in the original 
model, is determined internally at each geometric station. Also the friction 
coefficient relationship has been slightly changed. Comparison of results 
from the original and modified programs indicates that the performance degra- 
tion, the displacement thickness, and the temperature thickness are affected 
significantly. 

A study investigating the effect of important input parameters of the 
modified TBL version on the heat transfer and the performance degradation 
provides the following results: 

• A constant wall temperature or a wall temperature distribution 
shows a small effect on the heat transfer and performance degradation, if these 
temperature differences are moderate. 

• Consideration of a constant or internally calculated Prandtl number 
affects the heat transfer considerably since this quantity is contained in the 
Stanton number equation. The boundary layer thickness parameters are only 
slightly affected. 

• Use of two different Stanton number relationships produces a signif- 
icant heat transfer variation and affects boundary layer thickness parameters 
considerably. 

• Variation of the velocity profile exponent shows a definite effect on 
the boundary layer thicknesses and the performance degradation. The larger 
the exponent 1/n , the larger is the influence on performance. The heat 
transfer is hardly affected by the velocity profile. 


During this study a rocket thrust chamber using liquid oxygen and 

liquid hydrogen at a chamber pressure of P = 3000 psia ( 210. 9 kgf/cm 2 ) 

c 

and a mixture ratio of MR = 6. 0 , expanding the reaction products to an area 
ratio of e = 57 , was used as a reference case. 


INTRODUCTION 


In the original TBL program the boundary conditions of the boundary 
layer are represented by the wall temperature distribution along the chamber 
contour and by free stream flow properties correlated with a Mach number 
and geometry relationship along the boundary layer edge. Since the flow in 
the thrust chamber is chemically reactive during the expansion process, the 
edge conditions are not adequately determined from the ideal gas analogy. 
Therefore, the original program was modified to use profiles of Mach number, 
static pressure, static temperature, velocity, and mean molecular weight from 
an external source [JANNAF Two-Dimensional Kinetic (TDK) Computer 
Program] as boundary layer edge conditions. The specific heat ratio and the 
gas constant are internally calculated from a provided specific heat-temperature 
relationship and the mean molecular weight input. The original constant 
Prandtl number assumption was also replaced by an equation to determine this 
quantity at every local station internally. In addition the skin friction coefficient 
equation has been replaced by a modified Blasius relation to improve the 
solution accuracy and to avoid multiple solutions that were encountered at 
Reynolds numbers greater than R = 10 4 with the original formulation. 

% 

To identify solution differences, a rocket thrust chamber, originally 
considered for the Space Shuttle vehicle, was selected as a sample case. The 
results obtained with the original and modified TBL programs indicate significant 
differences and are presented in this document. 

A performance degradation analysis with the TBL program depends to 
a great extent on the input data. Therefore, a parametric study was conducted 
to identify the sensitivity of performance degradation and heat transfer with 
regard to certain input parameters. The critical input parameters must be 
determined to a much higher degree of accuracy than the ones showing only 
a minor effect. 

In this study the individual effects of the wall temperature distribution, 
Prandtl number, Stanton number, and the exponent of the velocity profile on 
the various boundary layer parameters, heat transfer, and performance 
degradation are determined with the modified TBL program. The results 
are documented and compared with the reference case. 
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The documentation of the original TBL program [1] was supplemented 
by an additional document [2] providing a detailed calculation sequence, a 
specific description of each subroutine, and derivations of the important 
equations used in the TBL program. 

The modified TBL computer program called TBL-I is available at 
S&E-ASTN-PPB, NASA, Marshall Space Flight Center. 


MODIFICATION OF THE JANNAF TBL COMPUTER PROGRAM 


The original TBL computer program [1] provided by the JANNAF 
Liquid Rocket Performance Committee internally computes the velocity, 
pressure, density, and temperature in the free stream of the combustion 
products from an ideal gas relationship. The reaction products flowing in 
the nozzle at a high temperature still react chemically even in the divergent 
section of the nozzle and are not accurately represented by the ideal gas 
analogy. Therefore, the TBL program was modified to accept tables of 
static pressure, static temperature, velocity, and mean molecular weight 
in addition to the thrust chamber contour and Mach number tables. These 
table values can be obtained from the JANNAF ODK or TDK programs [3, 4] 
or from another external source. 


Density (p« ) 

The density of the combustion products p is calculated from the 
perfect gas relation at each axial location, 

oo 

where P^ , T^ , and $1Z are the universal gas constant, static pressure, 

static temperature, and the mean molecular weight of the combustion products 
specified by the input tables. The subscript °° signifies the free stream 
condition. Equation (1) indicates that the density is affected by the mean 
molecular weight and that the specific gas constant represented by the ratio of 
the universal gas constant and the molecular weight is changing throughout the 
thrust chamber. 
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Specific Heat Ratio (?«,) 

The specific heat ratio y - can be represented by a speed of sound or 

Mach number relationship , 


M 2 !PT 

OO 


( 2 ) 


or by thermodynamic properties, 



where M , U , and 

oo 7 OO 

heat at constant pressure, 
chamber. 


(3) 

C denote the Mach number, velocity, and specific 

poo 

respectively, at each location within the thrust 


Combining equations ( 2) and ( 3) provides an expression for the 
specific heat C , 

poo 


C 

P 


OO 


u 2 



M 2 1pT 

OO < 


9TC 


(4) 


This equation could replace the existing specific heat-temperature 

input table in the program. However, since the C -T input table is equiv- 

P 

alent to the previous relationship and to avoid extensive reprogramming, the 
input table option was retained. 

Another option exists to calculate the mean molecular weight 3ft at 
each local station from an input table identifying the specific heat C at 

each axial station, 



M 2 T 

OO OO 


U 2 


(5) 
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The modified TBL program calculates the specific heat ratio y ' from 

equation ( 3) using the mean molecular weight table input. The specific heat at 
constant pressure C is provided by the table input relationship C =f(T ) 

poo P°® 

as in the original TBL program. 


Prandtl Number (P r ) 

The Prandtl number appears in the empirical, relation of the Stanton 
number and by definition represents a quantity based upon laminar flow. When 
the flow passes through the throat section, considerable changes occur in the 
Prandtl number. Therefore, the constant Prandtl number assumption has been 
replaced by an internal calculation of this parameter. Based upon the definition 
of the laminar Prandtl number, 


P 

r 



X 

oo 


( 6 ) 


where and X^ represent the viscosity and the thermal conductivity of 

the combustion products, the following equation was derived (Appendix A) and 
incorporated into the program: 



(7) 


This relation clearly shows the influence of the specific heat ratio y which 

in turn depends on the mean molecular weight and the specific heat 
according to equation ( 3) . A typical curve of Prandtl number variation for 
a nozzle flow field will be shown later in Figure 3. 


Skin Friction Coefficient (Cf) 

During operation of the original TBL program , multiple solutions of 
the skin friction coefficient were encountered for a given Reynolds number 
R (Appendix C). To avoid future calculation difficulties, Cole's skin 

e e 
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friction coefficient relationship was replaced by a modified Blasius correla- 
tion, 


where the Reynolds number uses the momentum thickness as the characteristic 
parameter, 



p U 6 

oo 00 


(9) 


According to experimental data [ 5] , the constant in the previous 
equation should be 0. 018 instead of 0. 025, but to remain conservative the 
constant value of 0, 025 was retained. In Figure 1 the curves representing 
the three discussed relationships are shown. The modified Blasius equation 
yields slightly higher values for the skin friction coefficient than Cole's 
relationship and is single valued in the higher Reynolds number range. Use 
of the constant 0. 018 shifts the modified Blasius equation results to lower 
values. 



0.3 0.4 0.5 .6 .7 . 8 . 9 1.0 2.0 4.0 6.0 8.0 10 20 


Figure 1. Skin-friction coefficient as a function of Reynolds number. 
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The skin friction coefficient C f ^ used in the Stanton number 
equation is replaced by a similar relationship , 



( 10 ) 


where the Reynolds number uses the energy thickness as the characteristic 
parameter, 



P U $ 




(ID 


Boundary Layer Performance Degradation 

In a rocket nozzle performance prediction, the losses associated with 
the existence of a boundary layer are significant. The original TBL program 
provides only the necessary parameters to determine the thrust degradation. 
For convenience, this equation has been included in this program. Since 
performance and performance losses are quite commonly expressed in terms 
of specific impulse, approximate equations for the thrust chamber vacuum 
and sea level thrust were incorporated into the program. This modification 
provides the capability to determine the specific impulse and impulse loss 
caused by the boundary layer. 

Two possibilities exist to determine the thrust degradation; these are 
dependent on the nozzle contour given. 

1. If the real nozzle contour is given, it can be postulated that the 
inviscid flow field will be almost the same as the one obtained with the real 
contour corrected by the displacement thickness, because the latter value is 
small compared to the contour radii. In this case the total thrust produced 
by a one -dimensional inviscid flow can be calculated, 


F = ir(r - 6* cos a} 2 ( P + p U 2 )1 

vacuum V / V 00 00 00 /J exit 


( 12 ) 
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and 


F i i = [^(r - 6* cos o; 2 ) (p - P + p U 2 jl . (13) 

sea level [_ \ / \ °° a ^<*> °° /I 

vXll 

The corresponding mass flow rate is 

m = L ,<r- 6* cos a) 2 . (14) 

L -I exit 

In these equations r represents the real nozzle radius. 

I f the postulation of negligible flow field differences does not hold true 
because of a sizable boundary layer displacement thickness, an iteration on 
the inviscid flow field contour is necessary. 

2. If the inviscid flow field contour is given, the thrust produced by 
a one-dimensional flow field is, 


F = |*7rr 2 (p + p U 2 jl 

vacuum I \ °° r *> °o yj 

and 


(15) 


F . . = [ir 2 (p -P +p U 2 ) . (16) 

sea level \ °° a / J ' 

The corresponding mass flow rate is, 

“ = I Poo u « ™ I • < 17 ) 

*- J exit 

In these equations r represents the inviscid flow contour radius. 

The specific impulse for the one-dimensional inviscid flow can now 
be calculated, 


I 

sp 

vacuum 


F 

vacuum 

m 


(18) 
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and 


sp 


sea level 


sea level 
m 


( 19 ) 


To calculate the thrust more accurately, the U m at the nozzle 

exit can be obtained from a given mass flow rate m . Depending on the real 
or inviscid flow contour, the mass flow density is 


P U 

* rV\ ft 


m 

7r(r - 6* cos a) 2 


(r = real contour) 


and 


( 20 ) 


p U = — r (r = inviscid contour) . (21) 

' OO oo 7rr* 

Multiplication of p x U m by the boundary layer edge quantity U m 

provides a more accurate result than the one obtained from the one-dimensional 
model. 


The thrust degradation caused by the boundary layer effects can be 
determined from the following equation, as documented in Reference 2; 


AF. 


B. L. 


[ 2?rr Poo u „ 20cosa ] ^ "ir 


P u 

r oo OO 


( 22 ) 


exit 


where the displacement thickness 6* and the momentum thickness 9 are 
the only boundary layer parameters. All other terms represent quantities 
of the boundary layer edge at the nozzle exit. 


The corresponding specific impulse loss is 




m 


(23) 
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The modified TBL program contains equations ( 12) , ( 13) , ( 14) , ( 18) , 
(19) , (22) , and (23) since the real contour of the engine, used in the param- 
eter study, was available. 


COMPARISON OF RESULTS FROM THE ORIGINAL AND 
MODIFIED TBL COMPUTER PROGRAMS 


A former Space Shuttle engine design was used as a sample case 
( Fig. 2) to show the difference in the boundary layer solutions obtained from 
the original and modified TBL programs. The original TBL program 
internally calculates the pressure , temperature T^ , and velocity 

in the free stream, whereas table values (Fig. 3 and Table 1) for these 
quantities and in addition the mean molecular weight 3ft along the boundary 
layer edge were used in the modified program. Because of difficulties 
encountered with the original skin friction coefficient calculation, the 
modified Blasius relationship was applied in both programs. A comparison 
of the results, shown in Table 2, indicates that the modified program calcu- 
lates a specific impulse loss approximately 50 percent lower, a total heat 
transfer approximately 5 percent higher, a boundary layer displacement 
thickness approximately 100 percent higher, and a momentum thickness 


X (ml 

0 12 3 



STARTING POINT 

Figure 2. Thrust chamber nozzle contour. 
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Xlm) 


0 1 I 3 



Figure 3. Inviscid nozzle flow properties as a function of 
nozzle length ( reference case) . 

approximately 25 percent smaller as compared to the original TBL solution. 

* All other solutions indicated only a slight deviation. The vacuum thrust 
calculation results differ significantly because of the internally calculated 
and tabular input nozzle exit parameters for the boundary layer edge. This 

• comparison clearly indicates that considerable differences result in the 
performance degradation, heat transfer, and boundary layer parameters 
when the modified program is used. 
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TABLE 1. INPUT TABLES FOR MODIFIED TBL PROGRAM 
(REFERENCE CASE) 


X 

(ft) 

y 

(ft) 

Mach 

Number 

Pressure 
( lbf/ft 2 ) 

Static 

Temperature 

r R) 

Velocity 
(ft/ sec) 

Molecular 

Weight 

(Ibm/lbm 

mole) 

Wall 

Temperature 

Cr) 

- 0.74000000 

1. 28300001 

0. 06540000 

430938. 589844 

6514, 171265 

341. 998299 

13. 596038 

1500. 000000 

- 0.70000000 

1.24500000 

0.08680000 

430803. 386719 

6513. 945374 

449.713486 

13. 596115 

1501. 000000 

- 0.65000000 

1. 15400000 

0.11500000 

430376. 687500 

6513. 232483 

601, 322884 

13. 596359 

1503. 000000 

- 0. 60000000 

1. 04400000 

0. 14800000 

429571. 242187 

6511. 884766 

773. 783501 

13. 596819 

1504. 000000 

- 0.55000000 

0. 93800000 

0. 18000000 

428257.367187 

6509. 682068 

940. 903488 

13. 597570 

1505. 000000 

- 0.50000000 

0. 84000000 

0. 21900000 

426143. 968750 

6506. 125793 

1144. 403244 

13. 598783 

1507. 000000 

- 0.45000000 

0. 74800000 

0. 26000000 

422584. 132812 

6500. 097595 

1357. 922546 

13. 600837 

1510. 000000 

- 0.40000000 

0. 67000000 

0. 30500000 

417121.234375 

6490. 756226 

1591. 619995 

13. 604016 

1512. 000000 

- 0.35000000 

0. 60000000 

0. 35500000 

408201. 140625 

6475. 257080 

1849.983658 

13.609287 

1515. 000000 

- 0.30000000 

0. 54000000 

0. 42000000 

394000. OOOOOO 

6449. 924011 

2183. 766510 

13.617876 

1517. OOOOOO 

- 0. 25000000 

0. 49300000 

0. 51000000 

373266. 746094 

6411. 375977 

2642. 587006 

13. 630889 

1522. 000000 

- 0.20000000 

0. 46200000 

0. 67500000 

348635. 917969 

6362. 908875 

3482. 375977 

13. 647145 

1525. 000000 

- 0. 15000000 

0. 44200000 

0. 83900000 

320962. 054687 

6304. 489929 

4305. 792297 

13. 666566 

1528. 000000 

- 0. 10000000 

0. 43000000 

0.91900000 

290185. 675781 

6233.701660 

4686. 352600 

13. 689814 j 

1530. 000000 

- 0.05000000 

0. 42500000 

0. 96900000 

257635. 841797 

6150.690002 1 

4904. 364929 

13.716626 

1535. 000000 

0. 00000000 

0. 42475000 

1. 00000000 

248237. 302734 

6124. 870544 

5049. 435852 

13.724856 . 

1540. 000000 

0. 26667000 

0. 60067000 

2. 02399999 

55255. 762695 

5103. 002380 

9304. 507935 

13. 986906 

j 

1532. 000000 

0. 70833000 

0. 94975000 

2. 70800000 

14754. 472656 

4211. 245239 

11400. 469604 

14. 087843 

1504. 000000 

1. 35833000 

1.20132999 

3. 03000000 

7766. 364746 

3790. 732147 

12161. 423950 

14. 103885 

1450. 000000 

1. 67500000 

1.34317000 

3. 18200001 

5748.760010 

3600. 483795 

12473. 540039 

14. 107473 

1415. 000000 

1. 95833001 

1.47142000 

3. 31400001 

4502. 654846 

3449. 974152 

12737.706665 j 

14. 109281 

1390. 000000 

2. 54166999 

1. 69900000 

3.51100001 

3068. 752136 

3221. 707977 

13073. 595337 

14. 110858 

1310. 000000 

3. 12500000 

1. 89950000 

3. 66999999 

2282.654144 

3052. 450104 

13327. 146729 

14.111449 

1243. 000000 

3. 91666999 

2. 12375000 

3. 83100000 

1699.312439 

2889. 791168 

13561. 844360 

14. 111748 

1163. 000000 

4. 66666996 

2. 32642001 

3. 96599999 

1336. 023575 

2761.740875 

13746. 858032 

14. 111872 

1093. 000000 

6. 58333004 

2. 68632999 

4. 18000001 

914. 612267 

2568. 120361 

14007. 482422 

14. 111960 

951. 000000 

8. 95833004 

3. 00342000 

4. 34299999 

681. 880005 

2424. 785431 

14171. 263184 

14. 111985 

820. 000000 

10. 19166994 

| 3.20675001 

4. 44400001 

573.919670 

2343. 315613 

14273. 130005 

14.111992 

780. 000000 
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TABLE 2. COMPARISON OF SOLUTIONS OBTAINED FROM THE 
ORIGINAL AND MODIFIED TBL PROGRAMS (REFERENCE CASE) 



Modified TBL 

Original TBL 

F , ton (lbf) 

207.658 

186. 603 

vacuum 

(457 810) 

(411 390) 

AF, ton (lbf) 

-1.587 

-2.072 

. 

(-3499) 

(-4567) 

Al , sec 

-3. 53 

-5.33 

sp 



q , kcal/cm 2 sec 

2. 802 

2.505 

w max 

( Btu/in. 2 sec) 

(71.7) 

(64. 1) 

Q , kcal/sec 

2. 2574 x 10 4 

2. 1505 x 10 4 

w 

(Btu/sec) 

(89 521) 

( 85 622) 

At nozzle exit 



6* , cm (ft) 

1.92 

0.88 


(0. 063) 

(0.029) 

6 , cm (ft) 

6.00 

6.98 


(0. 197) 

(0.229) 

A , cm (ft) 

9. 14 

9.33 


(0.300) 

(0.306) 

6 , cm (ft) 

0.46 

0.61 


(0.015) 

(0. 020) 

6*/0 

4. 150 

1.448 

<f> , cm (ft) 

0.76 

0.88 


(0.025) 

(0.029) 

P , T , U 

OO OO OO 

Input 

Internally 

calculated 




Sonic Point Start 


The original TBL program has the option of a sonic point start; i.e. , 
the computation of the boundary layer can be started at the nozzle throat. This 
option, however, is limited to the adiabatic wall case. Furthermore, the 
assumptions are made that the temperature and velocity thicknesses are equal 
and the derivative of the momentum thickness is zero at the throat location. 

Results from the modified TBL program indicate that the temperature 
and velocity thicknesses in the throat area are not equal (Figs. 4 and 5) . The 
solutions of the original TBL program utilizing the sonic start point option and 
the ones obtained from the modified program starting the calculation in the 
chamber show remarkable differences as indicated in Table 3. 

EFFECT OF SIGNIFICANT PARAMETERS ON 
PERFORMANCE AND HEAT TRANSFER 


In the following paragraphs, the effect of some significant parameters 
on the performance degradation and heat transfer is discussed. The modified 
TBL program was used for this investigation. Parameters such as the wall 
temperature distribution, the Prandtl number, the Stanton number, and the 
velocity profile exponent were changed one at a time, and the results were 
compared with the reference case. The program input data for the reference 
case were obtained from Reference 6 and represent one-dimensional equilib- 
rium solutions for a former Space Shuttle engine design (Figs. 2 and 3 and 
Tables 1, 4, and 5) . 


Effect of Wall Temperature Distribution 

The reference case uses a wall temperature distribution as shown in 
Figure 6 and Table i. This distribution was established using Reference 7 for 
the throat and Reference 8 for the divergent nozzle sections. Discontinuties at 
the end points were eliminated by plotting the individual temperature distribu- 
tions and selecting points to obtain a smooth transition. To investigate the 
wall temperature sensitivity on performance degradation and heat transfer, 
this temperature profile was replaced by a constant temperature of T^ = 1600° R 

(889°K) along the wall. The results of the TBL solutions are shown in Table 6 
together with the reference case. Surprisingly the wall temperature variation 
has only a small effect on the performance degradation and heat transfer 
(approximately 4 percent lower for the constant wall temperature case) . 
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Figure 4. Effect of wall temperature on the temperature thickness. 
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TABLE 3. COMPARISON OF SOLUTIONS OBTAINED FROM THE 
MODIFIED TBL PROGRAM (REFERENCE CASE) AND THE 
ORIGINAL TBL PROGRAM SONIC START OPTION 



Modified TBL 

Original Sonic Point Start 

F , ton (lbf) 

207. 658 

202.337 

vacuum 

(457 810) 

(446 076) 

AF, ton (lbf) 

-1.587 

-0. 978 


( -3499) 

(-2157) 

Al , sec 

-3.53 

-2.23 

sp 



q , kcal/cm 2 sec 

2.802 

0 

w max 

( Btu/in. 2 sec) 

(71.7) 

(0.0) 

Q w , kcal/sec 
(Btu/sec) 

2.2574 x 10 4 

0 

(89 521) 

(0) 

At nozzle exit 



, cm (ft) 

1.92 

3.17 


(0.063) 

(0.104) 

6 , cm (ft) 

6.00 

7.07 


(0.197) 

(0.232) 

A , cm (ft) 

9.14 

0.82 


(0.300) 

. (0.027) 

6 , cm (ft) 

0.46 

0.37 


(0. 015) 

(0. 012) 

8*/e 

4. 150 

8.592 

<f) , cm (ft) 

0.76 

0.024 


( 0. 025) 

(0.0008) 

P 

Variable value 


= 0. 90 = Constant 

r 



r 

P , T , U 

00 OO OO 

Input 

Internally calculated 
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TABLE 4. INPUT DATA FOR MODIFIED TBL PROGRAM (REFERENCE CASE) 



I 
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TABLE 5. C -T INPUT TABLE USED IN THE ORIGINAL 
P 

AND MODIFIED TBL PROGRAMS 


I 

Specific Heat 
(Btu/lbm °R) 

Temperature 

(°R) 

1 

6.3078549-01 

7. 0000000+02 

2 

6. 6436788-01 

1. 63800004-03 

3 

7.0660707-01 

2.3433153+03 

4 

7.2625109-01 

2. 56812004-03 

5 

7.4233390-01 

2.7617400+03 

6 

7.6597598-01 

3.05245004-03 

7 

7. 8021657-01 

3.22170004-03 

8 

8. 1702141-01 

3. 6004838+03 

9 

9. 1021608-01 

4.2112448+03 

10 

1. 1800579+00 

5. 1030000+03 

11 

1.6420093+00 

6. 16269204-03 

12 

1.6721517+00 

6.2363194+03 

13 

1. 7100503+00 

6.3319866+03 

14 

1. 7398426+00 

6.4101000+03 

15 

1.7576111+00 

6.4581300+03 

16 

1. 7674131+00 

6.4851350+03 

17 

1.7707793+00 

6.4944977+03 

18 

1.7735372+00 

6. 50220004-03 

19 

1.7781478+00 

6. 5151568+03 

20 

1.7799212+00 

6. 6000000+03 
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TABLE 6. COMPARISON OF SOLUTIONS OBTAINED FROM THE 
MODIFIED TBL PROGRAM USING DIFFERENT WALL TEMPERATURES 
(REFERENCE CASE USES A VARIABLE WALL 
TEMPERATURE DISTRIBUTION) 



433.3°K<T < 855. 6°K 

w 

/ g \ 

T - 888.9° K = Constant 
w 

/ g \ 


(780 o R<T w < 1540°Rj 

1 T^ = 1600° R = Constant J 

F , ton (lbf) 

vacuum. 

207.658 
(457 810) 

206. 778 
(455 867) 

AF , ton (lbf) 

-1. 587 

-1.516 


(-3499) 

(-3343) 

Al , sec 
sp 

-3.53 

-3.39 

q , 

w max 

kcal/cm 2 sec 
(Btu/in. 2 sec) 

2.802 

(71.7) 

2,779 

(71.1) 

• 

Q , kcal/sec 

2.2574 x 10 4 

2. 1660 x 10 4 

w 

( Btu/sec) 

(89 521) 

(85 895) 

At nozzle exit 



<5* , cm (ft) 

1. 92 

2. 13 


(0.063) 

(0.070) 

6 , cm (ft) 

6. 00 

6. 13 


(0. 197) 

(0.201) 

A , cm (ft) 

9. 14 

9. 88 


(0.300) 

(0.324) 

6 , cm (ft) 

0.46 

0.46 


(0.015) 

(0.015) 

6*/d 

4. 150 

4.674 

<p , cm (ft) 

0.76 

0.82 


(0.025) 

(0.027) 
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The boundary layer parameters, except the boundary layer velocity thickness 
and the momentum thickness, are higher by approximately 8 to 10 percent for 
the constant wall temperature profile. It should, however, be stated that the 
performance and heat transfer sensitivity is more pronounced if lower constant 
wall temperatures are considered. In Figure 4 and Figures 7 through 9, 
temperature thickness A , variation of the specific heat transfer rate , 

displacement thickness 6 * > and shape factor <5 r '/ 9 , are plotted as a function 
of the axial nozzle length. 



X|ft) 

Figure 7. Specific heat transfer rate. 
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Figure 8. Effect of wall temperature on boundary layer 
displacement thickness. 





X (m) 



Figure 9. Effect of wall temperature on the boundary layer 

shape factor. 

Effect of Prandtl Number 

In the reference case the Prandtl number is internally calculated 
according to equation ( 7) and covers a range from 0. 91 at stagnation condition 
to 0. 80 at the nozzle exit. This equation was replaced by constant values of 
Prandtl number P^ = 0. 51 and P r = 0. 91 during the sensitivity study. The 

results are presented in Table 7, and indicate that the performance degrada- 
tion is hardly affected. The heat transfer rate differs significantly for the 

two cases of P =0. 51 and P = 0. 91 , whereas a negligible difference 
r r 

appears between the P^ = 0. 91 and the reference case. This naturally holds 

true for the total heat transfer rate and indicates a 5 percent higher quantity 

for the P = 0. 51 case compared with P = 0. 91 solutions. Although the 
r r 

heat transfer is only slightly affected when the variable Prandtl number is 
replaced by P = 0. 91 , a considerable difference might result for high area 

ratio nozzles with extensive nozzle lengths. In Figure 10, the total heat trans- 
fer rate Q for the reference case and for P = 0. 51 is shown as a function 
w r 

of nozzle length. 
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TABLE 7. COMPARISON OF SOLUTIONS OBTAINED FROM THE 
MODIFIED TBL PROGRAM USING DIFFERENT PRANDTL 
NUMBERS (REFERENCE CASE USES INTERNALLY 
CALCULATED PRANDTL NUMBERS) 



P = 0. 51 = Constant 
r 

P = 0. 91 = Constant 
r 

0. 80 < P SO. 91 

r 1 

Variables 

F ton (lbf) 

207.737 

207.651 

207.658 

vacuum 

(457 981) 

(457 793) 

(457 807) 

AF , ton (lbf) 

-1.615 

-1.586 

-1.587 


(-3560) 

(-3497) 

( -3499) 

Al ( sec) 
sp 

-3.59 

-3.53 

-3.53 

q , kcal/cm 2 sec 

3.096 

2.798 

2.802 

w max 

(Btu/in. 2 sec) 

(79.2) 

(71.6) 

(71.7) 

Q , kcal/sec, 

2.375 x 10 4 

2. 2519 x to 4 

2.25474 x lo 4 

w 

(Btu/sec) 

(94 200) 

(89 305) 

(89 519) 

At nozzle exit 




6* , cm (ft) 

1.92 

1.92 

1.92 


(0.063) 

(0.063) 

(0.063) 

6 , cm (ft) 

6.04 

6. 00 

6.00 


(0.198) 

(0.197) 

(0.197) 

A , cm (ft) 

9.48 

9.11 

9. 14 


(0.311) 

(0.299) 

(0.300) 

e , cm (ft) 

0.46 

0.46 

0.46 


(0.015) 

(0.015) 

(0. 015) 

6*/e 

4.058 

4. 154 

4.150 

<p , cm (ft) 

0. 82 

0. 76 

0,76 


(0.027) 

(0.025) 

(0.025) 

Prandtl number 

Input 

Input 

Internally 

calculated 

P = 0. 91 at 
r 

the stagnation 
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Figure 10. Effect of Prandtl number on the total heat 
transfer rate. 



Effect of Stanton Number 


The original and modified TBL programs both use the following 
empirical relationship to determine the Stanton number: 



To examine the influence of Stanton number the previous equation was replaced 
by an equation from Reference 9 approximately representing the von Karman 
analogy, 



for this investigation only. 

A total heat transfer rate increase as much as 16 percent resulted with 
regard to the reference case, although the performance degradation was 
insignificant. 


Effect of Boundary Layer Velocity Profile 

The velocity profile across the boundary layer is commonly expressed 
by a power law, 



( 26 ) 



The recommended value for n in the exponent is n = 7 . This number has 
been used for the reference case although n is a program variable. Reference 
5 concludes that n changes with the Reynolds number and pressure gradient, 
and varies between n = 5 and n = 9 for a Mach number range from M = 3. 85 
to M = 4, 8 . It is shown in Reference 10 that n = 8 and 5 for Mach numbers 
M = 5 and 8, respectively. In studying the velocity profile effects on the 
boundary layer solutions, a value of n = 5 was selected in the power law. 

The results for the reference case and using n = 5 are presented in Table 8. 

A specific impulse loss increase of approximately 8 percent occurs when the 
exponent factor n is changed from n = 7 to n = 5, whereas a slight decrease 
in the total heat transfer rate is experienced. Figure 5 and, Figures 11 through 
13 show the velocity thickness <5 displacement thickness 6* , momentum 
thickness 9 , and shape factor <5 * /d as a function of the nozzle length for 
the reference case and the one using n = 5 . These parameters are affected 
considerably by the velocity power law, especially the velocity thickness 6 , 
and the effect increases with nozzle length. 


RESULTS OF THE PARAMETER STUDY AFFECTING 
PERFORMANCE AND HEAT TRANSFER 


For convenience the energy thickness cf> and the Reynolds number 
based upon the momentum thickness R g as a function of the nozzle length 

are provided for the reference case in Figures 14 and 15. Variation of the 
Reynolds number in the throat section is very large because of the influences 
of density, velocity, and momentum thickness. 

In Figure 16 the performance degradation in terms of specific impulse 
is presented as a function of nozzle length. Results obtained with the original 
and modified TBL programs as well as changes caused by the variation of the 
wall temperature and velocity profile with respect to the modified TBL version 
are exhibited. The effects of the Prandtl and Stanton numbers are not specif- 
ically shown, since the curves fall almost on top of the reference case. (See 
Table 7 and the section entitled "Effect of Stanton Number".) 

It is evident that the largest variation in the performance degradation 
occurs when the free-stream parameters calculated internally by a perfect 
gas analogy are replaced by more accurate input tables. The wall temperature 
and to a greater extent the velocity profile are significant parameters and must 
be accurately determined for a performance degradation analysis. 



TABLE 8. COMPARISON OF SOLUTIONS OBTAINED FROM THE 
MODIFIED TBL PROGRAM USING DIFFERENT FACTORS n IN 
THE VELOCITY PROFILE EXPONENT (n = 7, REFERENCE CASE) 



n = 7 

n - 5 

F , ton (lbf) 

vacuum 

207. 658 
(457 810) 

208. 643 
(459 980) 

AF , ton (lbf) 

-1.587 

-1.720 


(-3499) 

(-3792) 

Al , sec 
sp 

-3. 53 

1 

CO 

00 

o 

q , kcal/cm 2 sec 

w max 

( Btu/in. 2 sec) 

2. 802 
(71.7) 

2.802 

(71.7) 

Q , kcal/sec 

2. 2574 x 10 4 

2.2490 x 10 4 

w 

(Btu/sec) 

(89 521) 

(89 190) 

At nozzle exit 



6* , cm (ft) 

1. 92 

1.71 


(0.063) 

(0.056) 

6 , cm (ft) 

6. 00 

5. 09 


(0. 197) 

(0. 167) 

A , cm (ft) 

9. 14 

7.32 


(0.300) 

(0.240) 

6 , cm (ft) 

0.46 

0.49 


(0. 015) 

(0.016) 

6*/0 

4. 150 

3.488 

cf> , cm (ft) 

0. 76 

0.76 


(0. 025) 

(0. 025) 

u 

JL _ (i Y' 

u /'yV /5 

U 

QO 

A - w 

A '*/ 














X (ml 
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Figure 13. Effect of velocity profile power law on the 
boundary layer shape factor. 

As to the heat transfer, a sizable difference has been experienced when 
the original TBL solutions are compared with the modified TBL results. The 
Prandtl and Stanton numbers affect the heat transfer significantly, whereas the 
velocity profile produces only a small variation. The selected constant wall 
temperature indicated only a minor effect on the heat transfer. However, lower 
wall temperatures might have a much larger influence. Table 9 represents a 
quantitative estimate of the influence of the parameters investigated on the 
performance, heat transfer, and boundary layer parameters. 


CONCLUSIONS 


The original TBL program has been modified to accept boundary layer 
edge conditions from an external source. In addition, several equations within 
the program were modified, and new ones were introduced. A former Space 
Shuttle engine design was used to investigate how the boundary layer solutions 
were changed by the modifications. Comparison of the results indicated 
drastic differences in the performance degradation, heat transfer, and the 
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<t> (cm) 



X (ft) 

Figure 14. Energy thickness as a function of nozzle length 

( reference case) . 
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Figure 16. Specific impulse loss caused by boundary 
layer effects. 

magnitude of boundary layer parameters. During a parametric study investi- 
gating the effect of significant parameters on heat transfer and performance, 
results from the modified TBL program were used as a reference. A com- 
parison of the results indicated the strong influence of the velocity profile, and 
to a minor degree the wall temperature, on performance. The heat transfer 
rate was primarily influenced by the Prandtl and Stanton numbers. The 
influence of the wall temperature depends on the magnitude of the temperature. 
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TABLE 9. QUANTITATIVE ESTIMATE OF THE INFLUENCE 
OF INVESTIGATED PARAMETERS ON PERFORMANCE, 
HEAT TRANSFER, AND BOUNDARY LAYER PARAMETERS 



p 

r 

C H 


-=(*y /n 

u „ W 

T 

w 

g 


Solution Differences 
Between Original 
TBL and Modified 
TBL 

F 

1 

2 

2 

2 

1 


4 

vacuum 








AF 

1 

2 

4 

2 

2 


5 

AI 

1 

3 

4 

2 

2 


5 

sp 








q 

2 

5 

4 

1 

1 


4 

w max 








Q 

2 

5 

4 

2 

2 


3 

w 








6* 

1 

3 

4 

3 

2 


5 

6 

1 

3 

4 

4 

1 


4 

A 

1 

3 

3 

3 

2 


2 

e 

1 

3 

4 

2 

1 


3 

<P 

1 

2 

2 

1 

2 


3 


Notes: Influence; 1 — negligible, 2 — slight, 3 — considerable, 4 — large, 
and 5 — extreme. 


36 
















APPEND IX A 

DERIVATION OF PRANDTL NUMBER 


Let us denote the specific heat at constant volume as (cal/gr °C) , 
then the specific heat for one molecule with mass, m, is mC^ (cal/°C) . This 
is the increase in kinetic energy of a molecule with a temperature rise of 1°C. 

When a temperature gradient, d'p'dx, exists in a gas, the temperature 
difference of two molecules that are separated by a distance of 2 i, is equal to 
21 (dT/dx). Therefore, the energy difference between these molecules is: 

- dT 

mC x2i — — . (A-l) 

v dx 


Consider the mean -free -path length which is the average distance a 
molecule may travel before it collides with another molecule. The above two 
molecules exchange their energies by such a collision. Averaging the total 
molecules moving at random, we can assume that they move an amount of one- 
sixth in each direction of the Cartesian coordinates x, y, and z. By assuming 
there are N molecules per 1 cm 3 of a gas and their mean molecular velocity is 
v (cm /sec) , the number of molecules moving in one direction across the area 
of 1 cm 2 in 1 second is Nv/6. These transfer energies, as stated above, and 
their total heat transfer rate (cal/cm 2 sec) is: 


X 


dT 

dx 


(A-2) 


which is called the Fourier law. The thermal conductivity is X (cal/cm sec°C) 
and has a certain value depending upon the substance; therefore, we obtain: 


X 


dT 

dx 


mC x 21 
v 


dT Nv 
dx X 6 


(A -3) 
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MOLECULE WITH HIGH TEMPERATURE 



COLLISION 


MOLECULE WITH LOW TEMPERATURE 


so that 

X = Nmvi C = ~ pvl C 
3 v 3 v 


(A -4) 


where p = Nm is the density of the gas. 

Another relation for the thermal conductivity is derived in Appendix B, 

* - § ^ c p • (A - 5) 

where y is the ratio of the specific heat at constant pressure to that at 

constant volume, y = C /C . 

’ ’ p' v 

The viscosity of a gas is also caused by molecular collisions. By 
considering the momentum exchange between colliding molecules in a gas 
along with velocity gradients, the viscosity coefficient p (gr/cm sec) can 
be easily derived as in the case of equation (A-4) , 

P = ^ pvi" . (A -6) 
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I 


I 

\ 

I By dividing both terms of equations (A-5) and (A -6) respectively, we obtain: 


fiC 

P 

X 


9y - 5 


(A -7) 




i 


i 


t 

i 


The left side of equation (A-7) is the definition of the Prandtl number, P . 
Rearranging equation (A-7) results in 




Using the ideal gas relationship, 


C - C 
P v 


£ 

m 




finally yields 


P = 
r 


°P + 4 3H 


(A —8) 
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APPENDIX B 

DERIVATION OF THERMAL CONDUCTIVITY 


The internal energy of a molecule is composed of a transverse move- 
ment, rotation, and internal vibration. Therefore, the specific heat, C v , 
is composed of the components C^, C^, and C^: 


C = C , + C + C 


vt 


vr 


vo 


(B-l) 


Since the mean velocity of molecules is higher for those with higher 
temperature, the energy transfer by the transverse movement is large. 
According to Chapman 1 A is as much as 2. 5 times as great when molecules 
have only the degree of freedom of transverse movement instead of rotation 
or vibration only. Eucken 2 considered that influences of the rotation and 
internal vibration on the energy transfer are small even when both exist at 
the same time; then, equation (A-4) is 


A 



(2. 5 C + C + C ) 
vt vr vo 


(B-2) 


Substituting the following relations into equation (B-2) , we obtain 
equation (A-5) . Equation (B-l) becomes after rearrangement 


C +C = C - C . . 

vr vo v vt 


Using the gas dynamics relation, 



3 

2 911 ; 


(B-3) 


(B-4) 


1. Chapman, S.: Phil. Trans., vol. A211, 1912, p. 459. 

2. Eucken, A.: Forsch. Gebiete Ingenieurw. , vol. 11, no. 6, 1940. 
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and the ideal gas formulation, 



(B-5) 


In equation (B-2) , the following result for the thermal conductivity is 
obtained: 



1_ 9 y - 5 
3 4y 


pvi C p . 


(B-6) 
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APPENDIX C 

SKIN FRICTION COEFFICIENT 


The skin friction coefficient calculation in the original TBL program 
is based upon the following computation loop [ 11] . Initially an adiabatic 
skin friction coefficient C f is assumed. Then a corresponding skin friction 

a _ 

coefficient for low speed flow C ^ is calculated using Cole's relationship. 


P P c 


C f = 


P P 
r aw aw 


(C-l) 


where the subscript aw refers to adiabatic wall conditions, p is the free 

stream density, and ju represents the dynamic viscosity of the laminar 

s 

sublayer. The dynamic sublayer viscosity p is obtained using Cole's 

s 

temperature relationship. 


= 1+17.2 (|*- 
T \ T 

aw \ aw 



— I - 305 


aw 


aw 


2 ’ 
(C-2) 


in connection with 


P, 


m 


P 


aw 


aw 


(C-3) 


where the exponent n denotes the temperature dependence of viscosity. 
The following expression is then calculated according to Cole: 


R 


_PJL 


e e 


P M 
'aw ^aw 


R 


(C-4) 
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This equation is similar to the previous one except that the free stream 
dynamic viscosity M and the Reynolds number based upon the momentum 
thickness, 



pue 

M 


(C-5) 


are used. Now another low speed flow skin friction coefficient can be 
determined using 


C 


f 


0.009896 



(C-6) 


for a range of C R = 1 to 2. 51 , a tabular relationship based on test 
f 

data over the range of C f R = 2. 51 to 64. 8 , and 

1 e e 



for values of C. R > 64. 8 . 
f e » 

The resulting skin friction coefficient is then compared to the 

coefficient assumed initially. If both values of fall within a small 

tolerance, is satisfactorily determined, otherwise the calculation must 

a 

be repeated with an improved assumption. 

a 
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Application of the TBL program to a former Space Shuttle engine 
design, which was selected as a reference case in this document, observed 
that a plot representing the adiabatic skin friction coefficient as a function 
of the Reynolds number, showed multiple solutions for a constant Reynolds 
number for R g ^ > 10 4 (Fig. 1). Test data, however, indicate that such 

condictions do not exist. Therefore the previously described solution method 
for the skin friction coefficient was replaced by the following straightforward 
method. At first the Reynolds number at a specific location is calculated by 


R 

e 


0 


pue 

M 


(C-8) 


Then the skin friction coefficient is directly determined from the relation, 


c f ‘ 


K 


(\F 


(C-9) 


According to experiments by D. Brott [5] the constant K = 0. 018 for Mach 

numbers M > 1 with pressure gradients in the flow. In the modified TBL 

program a constant K = 0.025 was used instead representing more closely 

the Blasius relationship. From Figure 1 it is evident that solutions for the 

skin friction coefficient very closely compare with the original TBL concept; 

however, multiple solutions for Reynolds numbers R. > 10 4 are avoided. 

e 0 
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